Crossover from Isotropic to Directed Percolation 



\0 . 
0^ ■ 

a^ ■ 

Oh: 

^ ■ 

O • 
cn 



> 

oo 
(N 
0^ 
O 

ON 

a 

o 
o 



Per Frojdh and Marcel den Nijs 
Department of Physics, University of Washington, P.O. Box 351560, Seattle, Washington 98195-1560 

Directed percolation is one of the generic universality classes for dynamic processes. We study the 
crossover from isotropic to directed percolation by representing the combined problem as a random 
cluster model, with a parameter r controlling the spontaneous birth of new forest fires. We obtain 
the exact crossover exponent j/op = j/t — 1 at r = 1 using Coulomb gas methods in 2D. Isotropic 
percolation is stable, as is confirmed by numerical finite-size scaling results. For -D > 3, the stability 
seems to change. An intuitive argument, however, suggests that directed percolation at r = is 
unstable and that the scaling properties of forest fires at intermediate values of r are in the same 
universality class as isotropic percolation, not only in 2D, but in all dimensions. 



Directed percolation (DP) has emerged during recent 
years as one of the most common dynamic universality 
classes. It applies to a wide array of dynamic processes, 
ranging from flow through a porous medium in a gravi- 
tational field, forest fires and epidemic growth, to surface 
chemical reactions . In 1 + 1 dimensions the DP critical 
exponents are known accurately from numerical studies 
in the early eighties Today the numbers have been 
refined but analytic insight is still lacking. The ulti- 
mate goal is to understand scale invariancc in DP at the 
same level as that of isotropic percolation (IP). This will 
require some sort of generalization of conformal invari- 
ancc 0] and the Coulomb gas method |^,||, which apply 
to isotropic scaling phenomena in 1 + 1 dimensions. The 
scaling properties of dynamic processes like DP are in- 
trinsically anisotropic. Studies of crossover phenomena, 
like the one presented here, are a start into this direction. 

Consider a square lattice, as shown in Fig. 1. For IP 
the bond between vertices i and j represents a channel 
that can be open or closed, s.y = ±1, with a probability 
p or 1 — p, respectively. At p^— ^, an infinite cluster of 
connected channels appears, such that fiuid can perco- 
late all the way through the lattice from one end to the 
other. The correlation length diverges as ^ ~ \p — PcC^ 
with 1/v = yp = yr = 3/4 [||. 
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Directed percolation describes the same type of fiow 
in a porous medium, but in the presence of a gravita- 
tional field. The fluid can only flow downwards (posi- 
tive time direction in Fig. 1). Points at either end are 
considered connected only if there exists a path between 
them without the need to back-track in time. The DP 
threshold is larger, pc = 0.6447, and the scaling proper- 
ties at Pc are anisotropic. The time-like (||) and spatial 
(_L) correlation lengths diverge with different exponents, 
^11 ^ |p — Pel""" and £^±_ ~ \p — Pc\~''^ , respectively, with 
z = i'\\/v±_ ~ 1.581 the dynamic critical exponent and 
yp = l/j^j_ = 0.9117 the scaling dimension of p — pc- The 
values of these exponents are known only numerically . 

One of the alternate incarnations of DP is as a model 
for the propagation in time of a forest fire on a {D — \)- 
dimensional lattice. DP treats burnt trees as equivalent 
to never-burnt trees. This absence of immunization is 
unrealistic, unless we visualize this as a fire inside a D- 
dimensional forest in a strong wind. The fire propagates 
only in the direction of the wind, which maps out time. 
The difference between IP and DP is the possibility to 
spontaneously ignite new forest fires. Define two opera- 
tors (see Fig. 1): 
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equals one if a new forest fire is created at vertex i and 
is zero otherwise. 

A, = 4(1 - Sa){l - Sb){l " Sc)(l - Sd) (2) 
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equals one if vertex i is not part of any forest fire and is 
zero otherwise. Consider the following partition function 
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FIG. 1. Bond percolation on a square lattice with a 
time-like vertical direction: (a) Open bonds {sij = 1) are 
marked with solid lines, (b) Bonds around each vertex i are 
labeled Sa, Sb, Sc, and Sd- 



where iV^, is the number of vertices of the lattice, N^, 
the number of percolating bonds in graph Q, and the 
number of disconnected clusters in Q (including the ones 
that contain only one vertex). < r < 1 controls the 
spontaneous birth of forest fires, and the factors 
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1 — rp{2 — p) 



(4) 



ensure that the partition function remains stochastic at 
q = 1 for all values of r. At r = 1 we recover the conven- 
tional random cluster (RC) model |Q, which is equivalent 
to the g-state Potts model and describes IP in the q = 1 
limit 1). 

At r = 1 the RC model can be mapped onto the 6- 
vertex model and in the continuum limit onto a sine- 
Gordon type model. The latter can be rephrased as a 
Coulomb gas between charged plates [||. These equiva- 
lences lead in the early eighties to the derivation of the 
exact values of the critical exponents of the Potts model 
1^,^. The same techniques can be applied to determine 
the stability of IP with respect to r. At r = 1 and p — \ 
the crossover operator takes the form 

Oi = \{Sa + Sb + Sc + Sd) + \{Sc + Sd- Sa ~ Sb) 
o o 

- i [ScSd + 2{Sa + Sb){Sc + Sd)] + 0{s^). (5) 

First we need to substract from Eq. (5) operators that 
are already present in the theory at pc and r = 1 . Those 
operators are invariant under time-reversal T (sa <-> Sc 
and Sb Sd] see Fig. 1). For example, the spin variables 
Sij represent the energy density of the Potts model. It 
is invariant under T and changes sign under the duality 



transformation (s^ 



The first term on the right 



hand side of Eq. is the energy density. Its presence 
indicates that the critical line approaches r — I with a 
finite angle. The second term has the form of the gra- 
dient of the energy operator with respect to time, and 
breaks time reversal. This is the true crossover operator. 
It is essential not to miss any topological operators that 
might be hiding in Oi . This is where the mapping to the 
Coulomb gas becomes important. We have determined 
the Coulomb gas representation of the crossover-operator 
pair-correlation function. The details are not presented 
here, since we find no hidden topological operators. The 
most dominant component in Oi that breaks time rever- 
sal is simply the gradient of the energy operator with 
respect to time. 
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Hence, the crossover exponent is equal to i/dp ~ Vt 
This is an exact result. The crossover exponent is irrele- 
vant, hdp — —1/4 is negative since j/t — 3/4. This indi- 
cates that the scaling properties of a single forest fire, DP 
at r = 0, are different from those of forests in which new 
fires can spontaneously ignite, with probability propor- 
tional to r. Those must have the same scaling properties 
as IP. The following numerical results confirm this. 

We have performed a finite-size scaling (FES) analysis 
of Eq. (||) at q = 1. The transfer matrix acts on a state 



vector which in its most efficient representation contains 
only slightly less than 5^ elements for system size N . 
Non-local connectivity information is required at each 
moment in time. It does not suffice to know which trees 
are on fire and which are not. We need to know which 
burning trees at time t are part of the same cluster via 
a path in the past. Our coding for this is inspired by an 
earlier FSS analysis at r = 1 by Blote and Nightingale [|j . 
The full connectivity can be reconstructed by attributing 
to each vertex 5 possible states: (0, <^=, =>, ★). Vertex 
i is not part of any forest fire (state 0) if both trees just 
above it are not burning, Sa ~ Sb = —1 (see Fig. 1), but 
is part of the fire if either or both are burning. In state 
■i= at least one vertex to the left of i is part of the same 
cluster, but none to its right. In state at least one ver- 
tex to the right is part of the same cluster, but none to 
its left. In state <^ at least two vertices, one on each side, 
is part of the same cluster. In state * no other vertices 
at this moment in time are part of the same cluster. 
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FIG. 2. Finite-size estimates for the percolation threshold 
Pc{r) from the locations of the minima in the energy mass 
gap. Results for system sizes A'^ = 3, . . . , 8 (solid lines from 
bottom to top) are extrapolated to infinite A'^ (dashed line). 




FIG. 3. Finite-size estimates ^(A'^, A'^ -I- 1) for the dynamic 
exponent z using the relation m ~ A/N'. For r > 0, 
2;(3, 4), . . . , z{7, 8) (from top to bottom) converge to 2: = 1. 
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Fig. 2 shows FSS estimates for the location of the per- 
colation threshold pdr) for system sizes 3 < TV < 8. The 
leading and next-leading eigenvalues of the transfer ma- 
trix give us the energy mass gap, rrie = ^|~^, the inverse 
of the time-like correlation length of energy-energy cor- 
relations. This mass gap is finite at either side of Pc- The 
minima in rrie with respect to p at each r shown in Fig. 2 
converge to Pc{r). 

Fig. 3 shows FSS estimates for the dynamic exponent 
z. At Pc the energy mass gap scales with system size as 
nie — A/N^ . The values of the mass gap at the minima 
shown in Fig. 2, yield for each set of system sizes (N,N-Hl) 
a FSS estimate for z. These estimates converge to the IP 
value z — 1 for all r, except very close to r = 0. Crossover 
scaling toward the DP value z = 1.581 does not set in 
until r is rather small. FSS estimates for yp = l/v^ (not 
shown here) can be obtained from the derivatives of rrie 
with respect to p at the minima. These behave similarly 
to the estimate for z in Fig. 3 and converge towards the 
IP value i/p = j/T = 3/4 for all r > 0. Crossover scaling 
towards the DP value does not play a role until rather 
small values of r. These numerical results clearly confirm 
that for all ?■ > the scaling properties are identical to 
those of IP. 
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FIG. 4. Finite-size estimates for the anisotropy parameter 
fj, = C-l/Cii s-t Pc using the relation m ~ 2nfj,{2 — yT)/N for 
system sizes — 3, . . . , 8 (solid lines from top to bottom). 
The dashed line corresponds to the extrapolation N ^ oo. 

Fig. 4 shows the FSS estimates for the aspect ratio 
/i = at Pc{r). Both correlation lengths diverge 

with exponent v — 4/3, but their amplitudes differ by 
a factor /x. We obtain from the energy mass gap. It 
scales according to conformal field theory as nie — A/N 
with a universal amplitude A = 27r/Li(2 — yx)- In Fig- 4 
we plot Nme/2TT{2 — yr) (with yx = 3/4) at the min- 
ima shown in Fig. 2. For DP this ratio represents the 
so-called opening ang le § and goes to zero at Pc- The 
convergence to a non-zero value for r > gives further 
confirmation that the scaling is isotropic, z = 1. 

It is a safe guess that Eq. (P) and yop ^ yr ~ hold 



also in dimensions larger than D = 2. Eq. is easily 
generalized to higher dimensions. The only aspect miss- 
ing for general D is the Coulomb-gas rigorous proof that 
the scaling dimension of the gradient of the energy main- 
tains its naive value yr — ^- The relation yop = Vt — ^ 
suggests that IP changes stability with dimensionality. 
The critical exponents for IP in 3D are known numeri- 
cally 0; yr = 1.14 yields yop = 0.14. So for D > 3 
the crossover operator appears to be relevant. Mean- field 
theory, valid in D > 6, predicts that the crossover expo- 
nent is relevant as well. Power counting yields the value 
VDP = 1. 

Our identification of the crossover scaling field with the 
energy gradient agrees qualitatively with a recent field- 
theoretical analysis by Frey et al. Q (see also Ref. @). 
They describe IP by a (/)^-type theory similar to the con- 
ventional one for the g-state Potts model, and set up 
an interpolation scheme between this and Reggeon field 
theory which belongs to the DP universality class. 
Frey et al. perform a one-loop expansion about the up- 
per critical dimension = 6 for IP and observe indeed 
a flow from isotropic to directed percolation in their field 
theory. Our approach has the advantage of being based 
on a well-defined microscopic model, and being exact in 
2D. 

The apparent turn-about in the stability of IP for 
2 < 13 < 3 is the most intriguing implication of Eq. (|^). 
Initially wc thought that IP should be unstable towards 
DP in all dimensions, r < 1 introduces a preferred direc- 
tion in the model (the time- like direction). This breaks 
a fundamental symmetry of the problem and therefore 
it seems logical that the scaling properties change. This 
is too simplistic, however. The crossover operator is the 
gradient of the energy density, j/t is the fractal dimension 
of so-called "red vertices" [Q, i.e., of vertices where the 
cluster is singly connected (like the * states in our trans- 
fer matrix). Such vertices are rare, to the extent that in 
D — 2 their fractal dimension becomes smaller than one, 
yr — 3/4. This causes yup to become negative. 

The following argument puts the shoe on the other 
foot. It supports the opposite point of view, namely that 
DP at intermediate values of r belongs to the same uni- 
versality class as IP (at r = 1) in all dimensions. Consider 
DP at r = just below its percolation threshold. All for- 
est fires that are initially burning die out. They are finite 
in size with a characteristic aspect ratio /io = At 
small non-zero r, new forest fires are starting at all times 
t. They have the same typical size and aspect ratio as 
those at r = 0. The probability that a specific vertex be- 
longs to a forest fire is of order p ~ rt^~^i\y These forest 
fires are placed at random, and increasing r is equivalent 
to increasing p in a conventional percolation problem, one 
in which objects of a specific size and shape are placed 
at random in space. Let pc be the percolation threshold 
for that percolation problem. The scaling properties at 
Pc are identical to those at r = 1. 
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The above argument suggests that DP is always un- 
stable with respect to the creation of spontaneous for- 
est fires, ft is probably too simplistic to treat pc as a 
constant, but doing so leads to the following estimate 
for the crossover exponent yip at r = 0. The IP crit- 
ical line approaches the DP critical point at r = as 
^ ^^±"^^11 r{p - pJ-(^+^-i)''^ with Pc the DP 
percolation threshold at r = 0. According to conven- 
tional scaling theory the critical line must approach the 
DP critical point along a renormalization flow line. This 
yields yip = z + D — 1 since yp = l/v±. This value can 
only be an approximation, since it violates the funda- 
mental rule that critical dimensions, like y/p, cannot be 
larger than the embedded dimension. It shows however 
that DP is strongly unstable with respect to r for all D. 

This leaves us with a puzzle. Do the scaling properties 
at intermediate values of r change between D — 2 and 
D = 3, as suggested by the stability analysis at r = 1, 
Eq. (|^), or not, as suggested by the above argument close 
to r = 0? Future numerical studies in D — 3 will decide 
this issue. For the time being, we put our bets on the 
above argument. It is hard to see how it can be circum- 
vented. Moreover, the gradient nature of the crossover 
operator at r = 1 is suspicious. Gradient operators can 
be integrated up into boundary operators and then often 
vanish from the theory altogether. This would leave the 
crossover at r = 1 irrelevant in all D. 

This research is supported by NSF grant DMR- 
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